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ABSTRACT 

We are presenting a wide-field image of the Mon R2 star forming region 
obtained with WFCAM on UKIRT in the 2.12 /im filter centered on the H 2 1-0 
S(l) emission line. We report the discovery of 15 new H 2 jets in Mon R2 and 
two in L 1646 and confirm most of these discoveries using archival Spitzer IRAC 
4.5 /im and 8.0 /im images. We find that many of these protostellar jets are found 
in projection against the outflow cavities of the huge CO outflow in Mon R2, 
suggesting that the jets may be associated with an episode of star formation 
in Mon R2 triggered by this large, but now fossil, outflow. We also study the 
spatial distribution of small, localized reflection nebulae and find that these are 
distributed in the same way as photometrically identified Class I sources. 

Subject headings: stars: pre-main-sequence — stars: formation — ISM: jets and 
outflows — ISM: Herbig-Haro objects 



INTRODUCTION 



The region studied here was first identified by Ivan den Bergh I (119661 ) as an association 
of optical reflection nebulae in their list of such objects and was labeled Monoceros R2, 
in short Mon R2. The Mon R2 molecular cloud complex extend s over 3° x 6°, and the 
overall mass of this large molecular cloud has been estimated by iMaddalena et al.l (119861 ) 
to be 9 x 10 M . The Mon R2 molecular cloud contains several sites of star formation: 



GGD 11, GGD 12-15 (IGvulbudaghian. Glushkov. fc Denisvuk Ill978f ). L 1646 flCarpenter 



2000|), and the Mon R2 core, the main star formation site that is the subject of this study. 
The Mo n R2 star forming reg ion lies between the two optical reflection nebulae vdB 67 and 
vdB 69 (Ivan den Bergh Ill966l ) and has attracted considerable attention as a site of massive 
star formation, as a young embedded cluster, and as the source of one of the most powerful 
CO outflows known. 
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Even though Mon R2 lies close to Or ion in projection on th e sky, its distance of 
830±50 pc based on photome tric parallaxes (Herbst &: Racine 1976h is almost double that 
of the Trapezium cluster, e.g. (jStanke. McCaughrean. fc Zinnecker II2002I ). Despite its larger 
distance, it is still close enough for detailed studies of the emb edded young object s. The mass 
of molecular gas in the Mon R2 core has been determined by lRidge et al.l (120031 ) using 13 CO 
and C 18 maps. The rarer C 18 traces denser gas and resulted in a gas mass of 1826 M , 
while the more abundant 13 CO also samples the tenuous outer regions of the molecular core 
and gave a mass of 2550 M Q . The molecular gas distribution is single-peaked with a slight 
indication of extended filaments to the south and south-west. To put this region into per- 
spective, Mon R2 is smaller than the Orion Trapezium region in spatial extent, molecular 
mass, and number and total mass of embedded stars. It is, however, not far behind in many 
of these parameters. 



Mon R2 has been noted as a source of a bipolar CO ou tflow by lLoren I ( 1l98ll ). This main 
Mon R2 molecular outflow was found by I Wolf et al.l (119901 ) to be one of the dynamically old- 
est (1.5xl0 5 y), largest (7 pc), and most massive (1 00 M P) ) molecular o utflows. Based on CO 
maps with higher spatial resolution than those of IWolf et al.l ( 1l990l ). iMeyers-Rice fc Lada 



( 119911 ) argued that the Mon R2 CO outflow must consist of two bipolar outflows: one compo- 
nent seen almost head-on to explain the overlapping red and blue shifted velocity components 
seen in projection against the central cluster, and one component more inclined against the 
line of sight, to explain the bipolar outflow features. Both components could be modeled 
as paraboloidal shells with a linear velocity field that could result from velocity-sorting of 
clu mpy outflow materia l. The detailed study of CO emission and multiple transitions of CS 
by iTafalla et al.l (119971 ) showed that the massive outflow has created an hourglass-shaped 
cavity in the molecular cloud centered on the cluster of massive young stars. The cavity 
walls coincide with the limb-brightened shells of the blue-shifted outflow lobe. The axis of 
the paraboloidal outflow shell and cavity are oriented roughly at a position angle of -30°. 

The central gro up of luminous stars a nd the associated nebulosity were discovered in 
the near infrared by lBeckwith et al.l (119761 ). Much of th e attention has b een fo cused on the 
ring-like reflection nebula around IRS 1 discovered by iBeckwith et al.l (119761 ) and studied 



in more detail by iHodappI (119871 ). on the details of the cluster's most 



source (IRS 3) summarized in the most recent work by iPreibisch et a 



uminous individual 
(l2002h. and on the 



distribution of molecular matter and the large outflow, summarized in ITafalla et al.l (119971 ) . 
In embedded young clusters like Mon R2, near-infrared (J, H, and K) photometry is a very 
useful tool for identifying young cluster members on the basis of their color excess over the 
re ddened photospheres of old background stars. The first such study of Mon R2 was done 
by ICarpenter et al.l (119971 ) . Most rec ently, the Initial Mass Function (IMF) in the Mon R2 



cluster core was studied in detail by lAndersen et al.l (120061 ) on the basis of HST/NICMOS 
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data. They find that the ratio of the number of stars below solar mass to the number of 
objects in the brown dwarf mass range is similar to that found in the Trapezium cluster and 
in IC 348. 

Compared to the widely used J, H, and K images used to characterize the stellar pop- 
ulation of young clusters, the Spitzer IRAC images in four bands ranging from 3.6 /im to 
8.0 /im allow the identification of younger and more deeply embedded objects. Further, 



( Gutermuth 


2005) and ( 


Kumar et al. 


2007) 



the basis of n ear- in frared J, H, and K photometry, and of Spitzer mid-infrared photometry, 



Gutermuth I (120051 ) found that Class I objects in Mon R2 are distributed in a filamentary 
pattern, while Class II objects are distributed more widely and evenly. However, even these 
IRAC bands are not sampling the youngest objects in a cluster completely, in particular 
those Class and I objects that are oriented with their disks seen edge-on. Fortunately, 
those objects can still be found indirectly. The substantial accretion activity in the Class 
and I phases of low-mass star formation is invariably associ ated with outflow activity, as 



was first noted by I Andre. Ward-Thompson fc Barsony I (119931 ). These outflows easily break 



out of the dense molecular disk and the envelop enshrouding the young star, and are there- 
fore readily detected in the near-infrared, due to the shock-excited emission at the interface 
between the outflow and the ambient molecular material, and in shocks within the outflow. 
The most easily accessible near-infrared shock-excited emission lines are the [Fell] line at 
1.644 jum and the H 2 1-0 S(l) line at 2.122 /zm, the line studied here. Of these two emission 
lines, the H 2 1-0 S(l) line traces lower excitation levels and is excited in virtually all deeply 
embedded jets propagating in molecular clouds. Therefore, shock-excited emission of molec- 
ular hydrogen is the earliest near-infrared signpost of outflow activity in a very young object 
and is usually observable long before the central star becomes observable in the near or mid 
infrared. Molecular hydrogen jets tend to have morphological features that are distinct from 
other forms of H 2 emission in star forming regions, and are therefore easily recognized in 
imaging surveys. 

Powerful jet activity is a byproduct of the short-lived main accretion phase of a forming 
star. When the accretion and jet activity fades after about 10 5 y, the envelope of molecular 
gas around a young star usually has formed two outflow cavities. When properly oriented, 
scattered light from these cavity walls gives the object the appearance of a bipolar nebula, 
or a cometary nebula if only one cavity can be seen. Objects of this morphology overlap 
with the Class I objects that can be identified photometrically using the Spitzer IRAC 
bands. Therefore, imaging surveys of star forming regions for protostellar jets and bipolar 
or cometary nebulae are a powerful tool for discovering very young (Class or I) objects in 
or soon after their main accretion phase. Thanks to the recent availability of large-format 
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near-infrared cameras such surveys can now be conducted efficiently. 

In this paper, we report the discovery of 17 extremely young objects in and around the 
Mon R2 star forming region by means of imaging their collimated outflows in the H 2 1-0 
S(l) emission line at 2.12 ^m. By comparison of the 2.12 /im image and the Spitzer IRAC 
4.5 /im and 8.0 fim images we discuss the most plausible identification of the source of each 
outflow. The large number of outflows indicates that star formation is active all through the 
Mon R2 molecular cloud. We will discuss the spatial distribution of the newly discovered 
extremely young stars in relation to the large, old molecular outflow in Mon R2. 



2. OBSERVATIONS AND DATA REDUCTION 

2.1. UKIRT WFCAM 2.12 /an Imaging 

The purpose of this project was not the detailed study of any particular jet, but rather 
the discovery of new objects in the Mon R2 star forming region that indicate the presence 
of stars in their youngest phase of evolution. Also, completeness in the discovery of H 2 
outflows was not a realistic goal for this survey of a region with such extensive reflection 



nebulo sity and, most likely, also fluorescently excited H 2 emission (e. g.. lBlack & Dalgarno 



( 119761 )) that both tend to confuse any search for H 2 outflows. Therefore, in the interest of 
maximizing the scientific return of the observing time used, we did not obtain an image at 
a continuum wavelength adjacent to the S(l) line, or in a broad K-band filter to distinguish 
between emission line features and continuum reflection nebulosity. As discussed below, the 
identification of H 2 jets relied on morphological arguments and confirmation, if possible, by 
archival Spitzer IRAC images. 

Wide-field near-infrared images of the Mon R2 region were obtained in the night of 
October 20, 2005, UT, at the United Kingd om Infrared Telescope (UKIRT). The Wide- 



Field Camera (WFCAM) flCasali et all 120071 ) was used with a filter centered at 2.121 /im 



with AA = 0.021 /mi that includes the H 2 1-0 S(l) line. The field was roughly centered 
on the well-studied Mon R2 cluster of infrared sources. The WFCAM uses four Teledyne 
(formerly Rockwell) HAWAII-2 2048x2048 HgCdTe infrared detector arrays at a pixel scale 
of 0.4" pixel -1 , arranged in the camera focal plane with a spacing of 0.94 of the detector size 
between the devices. To achieve full coverage of a 0.75 square degrees region of the sky, four 
telescope pointings are required, which in WFCAM terminology are called a "tile" . At each 
tile position, a five-point dither pattern (3.2" spacing) is used to avoid bad pixels and other 
detector array artifacts. Finally, at each dither position, a four-point micro-step pattern 
with a half-pixel (0.2") step size is used to oversample the image. For our observations, an 
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on-chip exposure time of 40 s was chosen, and the detectors were read out in non-destructive 
mode, sampling the accumulating signal 41 times during the integration time. The readnoise 
achieved in this mode of operation is ~ 20 e~ rms. The on-sky integration time was 800 s, 
except in the overlap regions of the tile. With a distance to Mon R2 of 830 pc, one arcminute 
in our images corresponds to a projected distance of 0.24 pc. 

The data were proce ssed by the Camb ridge Astronomical Survey Unit (CASU) using the 
procedu res described by iDve et al.l (120061 ). Astrometric calibration relative to the 2MASS 
catalog (ISkrutskie et al.l 120061 ) was performed and a plate solution was written into the file 
headers. All coordinates in this paper therefore are given in J2000.0 and are indirectly based 
on the 2MASS coordinate system. After retrieval of the reduced data from the WFCAM 
Science Archive the final assembly of the full tile image was done with IRAF scripts that used 
the image header information to assemble the individual images into a tangential projection 
image of the full tile. 

The central region of the resulting image is shown in Fig. 1. The image shows the well- 
known Mon R2 cluster near the center, and the two reflection nebulae to the west (vdB 67) 
and east (vdB 69) of it. Filamentary emission extends mostly to the north and east of 
th e main star formi ng region. The lowest contours of the blue and red shifted CO maps 
of IWolf et al.l (119901 ) are superposed on Fig. 1. There is no large-scale feature visible that 
would be morphologically associated with the main CO outflow centered on IRS 2 (identi fied 
in Fig. 2) that has a rather poorly defined outflow axis at P.A. ^-45° ([Wolf et al.lll990l) or 
with the outflow cavity whose axis is oriented at ~-30° (ITafalla et al.lll997l ) and (jChoi et al. 

2oooh . 



2.2. Archival Spitzer Images 



The IRAC camera is the main n ear and mid-infra red imaging instrument of the Spitzer 
cryogenically cooled space telescope (IFazio et al.ll2004j ). All the four IRAC channels contain 
emission lines of H 2 . Of the four IRAC channels, band 2 centered at 4.5 /im is particularly well 



suited to detect H 2 emission. From a discussion of different shock models, ISmith fc Rosen 



(120051 ) concluded that shock-excited H2 emission in the 4.5 /im band is an order of magnitude 
brighter than in the other channels for a wide range of shock conditions. The IRAC band 2 
(4.5 /im) contains the H 2 0-0 S(9) line at 4.694 /im and emission in this line shows morpho- 
logical features closely identical to those seen in H2 1-0 S(l) at 2.122 /im, but is less affected 
by dust extinction. Therefore, despite the fact that the broad IRAC band 2 also contains 
other important molecular features such as CO (v=l-0 at 4.45 - 4.95 /im) and the atomic 
Hydrogen Bra line at 4.052 /im, images in this band often show very similar jet morphology 
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to 2.12 /tm S(l) images (jSmith et al.ll2006l ). By contrast, ISmith &: Rosen I (120051 ) find that 
IRAC band 1 (3.6 //m) is dominated by vibrationally excited lines with higher excitation 
energy than the 1-0 S(l) line, so that shock features appear sharper in band 1 than in a 
2.12/zm 1-0 S(l) image. The Spitzer 3.6 /tm, 5.8 /tm, and 8.0 /tm bands also contain bright 
PAH emission at 3.3 /tm, 6.2 /tm, and 7.7 /tm, respectively. In a region like Mon R2, this 
strong and often filamentary emission tends to mask shock-excited emission from jets. 

Spitzer IRAC images were downloaded from the Spitzer archive. The data used here 
were originally obtained for a program by G. Fazio. Partial results of their study of the distri 



bution of stars based on JHK and Spitzer/IRAC colors have been presented by iGutermuth 
(£3). 



The available archival Spitzer data do not cover the whole field of our WFCAM ob- 
servations, so that we do not have the corresponding longer-wavelength data for all of the 
outflows discovered here. For those objects where all Spitzer bands are available, we discuss 
the 4.5 jttm and 8.0 /tm data. The 4.5 /mi data usually show the outflows, but at relatively low 
signal-to-noise ratio and with poorer spatial resolution that the UKIRT data. The 8.0 /mi 
data do not usually show the outflow itself, but in many cases show the molecular clump 
around the driving source either in thermal dust emission or in some cases in absorption 
against PAH background emission. The 8.0 /tm images therefore primarily serve to confirm 
the identification of the driving source. In some cases the band 2 (4.5 /tm) data were not 
available and for one object band 1 data (3.6 /tm) data are shown instead for the purpose of 
source confirmation. 



2.3. Criteria for Identifying H 2 Jets 



This simple imaging project was motivated by the succ ess of similar obs e rvatio ns in 
identifying H 2 jets, e. g. the e arly survey of NGC 1333 by iHodapp fc Ladd I (119951 ). the 
extensive survey of Orion A by IStanke. McCaughrean. &: Zinnecker I (2002). the survey of 
the distant massive star forming reg i on W 51 by IHodapp &: Davis I ( 120021 ) . and the recent 



work on DR21/W75 by iDavis et al.l (120071 ). Well over one hundred of these collimated H 2 
jets have now been imaged and their morphology is quite unique and well recognizable. Em- 
pirically, shock excited H 2 emission can, in most cases, be distinguished from fluorescently 
excited H 2 emission on the basis of morphological arguments. The physical basis for this 
i s that the cooling time of shock-excited molecular hydrogen is typically only a few years 
( jShull &: Hollenback Ill978l ). during which time shock fronts travel only a fraction of an arc- 
second in typical star forming regions. Shock-excited emission therefore usually appears 
as small but resolved knots, filaments, and bow shocks. In contrast, fluorescently excited 
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molecular hydrogen emission varies on the scale of typical density inhomogeneities in molec- 
ular clouds and therefore tends to be much more smoothly distributed. The same is true for 
large-scale continuum reflection nebulosity. 

Our criteria for identifying H2 jets in Mon R2 therefore were that they must consist of 
a string of resolved H2 emission knots. If a bow shock is found at the end of a string of 
emission knots, this is regarded as additional and convincing evidence for the existence of a 
jet. In cases where both sides (lobes) of a jet are visible, a high degree of central symmetry 
in the distribution of emission knots serves as supporting evidence for the identification of 
a jet. Sometimes, but not always, faint filamentary emission arising from the outflow cavity 
walls connects the individual emission knots. If such filamentary emission is seen, it serves 
as additional evidence supporting the identification of a jet. The H 2 jets identified by these 
criteria are listed in Table 1, and identified in Fig. 1 as large circles. 

Finally, many young stars show reflection nebulae of bipolar or cometary morphology 
or nebulae otherwise clearly associated with individual stars, properties that we summarize 
by the term "localized". In addition to their characteristic morphology, they show a more 
pronounced drop in surface brightness away from the illuminating central star than the shock 
emission knots in H2 jets. In Table 2, we list such localized nebulae that were found on our 
2.12/mi image on Mon R2, but were not classified as likely H2 jets. Nebulous features that are 
simply knots in the large-scale nebulosity permeating the center of the Mon R2 cluster are 
not listed in Table 2. Similar to the situation with the H 2 jets, this list cannot be expected to 
be complete, due to sensitivity limits, high extinction in some parts of the molecular cloud, 
and confusion with the large-scale reflection nebulosity and filamentary emission present in 
the Mon R2 region. The localized reflection nebulae are indicated as small squares in Figs. 1 
and 2. 



The K' survey of known CO outflow sources by iHodapp I (119941 ) showed that many of 
these outflow sources are associated with such localized reflection nebulae. Reflection nebu- 
lae, in particular of bipolar or cometary morphology, are produced by light scattered from 
the inside walls of the cavity produced by an outflow in the molecular material surrounding 
a young star. While the association of reflection nebulae with young age of an individual 
star and with outflow activity is less direct than for H 2 jets, the presence of large numbers 
of these localized reflection nebulae nevertheless indicates a population of very young stars, 
probably stars just after their main accretion and outflow phase. 
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RESULTS AND DISCUSSION 



3.1. The Environment of the Mon R2 Cluster 



Fig. 1 shows an extended system of emission filaments to the north and east of the 
cluster, besides the known reflection nebulae vdB 67 and 69, the nebulosity associated with 
the Mon R2 cluster, and the individual H 2 jets reported here. 

These filaments coincide with the optical reflection nebula vdB 68, and are also strongly 
visible in all Spitzer IRAC bands, probably indicating a combinatio n of H_2 and PAH emissio n 
fluorescently excited by the illuminating star of vdB 68, a B2V star (IHerbst fc Racine Ill976l ). 
The filaments seen in Fig. 1 are oriented at a position angle of ~ 35° , similar to some of 
the filamentary structure seen in the optical reflection nebula. It should be noted that this 
orientation of the filaments is similar to t he polarization angl e of background star polarization 
in the region of vdB 68 as measured by iJarrett et al.l (119941 ) in the R band, which indicates 
the projected direction of the local magnetic field. 

On the Digital Sky Survey (DSS) blue and red images, an extensive system of filaments 
is seen in absorption against the background reflection nebulosity in the area south and east 
of vdB 67, the Mon R2 cluster, and vdB 69. These filaments are also orien ted along the 



direct ion of the local magnetic field, as measured by the R-band polarimetry of IJarrett et al. 



(1l994f ) . Our 2.12 fim S(l) image does not show any emission from these filaments located in 
the south east of the embedded cluster. 

From the fact that these systems of filaments are readily visible at optical wavelengths, 
we conclude that they lie on the front side of the Mon R2 molecular core. Sinc e the polariza - 
tion vectors of more deeply embedded sources in and near the Mon R2 cluster (lHodapplll9871 ) 
follow the pattern of the magnetic field measured to the south and west of the cluster, and 
since the cluster is located at the intersection of the two patterns of magnetic field direction, 
it can be speculated that the formation of the Mon R2 clus ter was triggered by a collision 
to two molecular clouds, as discussed by IJarrett et al.l (jl994|). 



The core of the Mon R2 cluster and the bright reflection nebulosity permeating it are 
shown in more detail in Fig. 2. From polarization measure ments of this r ing-like infrared 
reflection nebulosity near the core of the Mon R2 cluster by iHodappi (119871 ). it is clear that 
IRS 2 is its dominant source of illumination. The brighter source IRS 3 does not significantly 
contribute to the illumination (as measured by polarization) of the ring-like nebula, and 
therefore must be located either in front or behind the nebula and the core of the Mon R2 
cluster. Some elongated emission features are seen east of IRS 3 and lie in a direction similar 
to the filaments observed near vdB 69. This would suggest that IRS 3 is illuminating some 
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of these filaments and that it lies in front of the rest of the Mon R2 cluster. 



3.2. Previously Known Outflows 



The overwhelming size and total mass of the main Mon R2 CO outflow (jLoren Ill98ll ) 
makes it difficult to find evidence for other o utflows in the velocity field m easured in 12 CO(2- 



1) or 13 CO(2-l). It should be noted that iMeyers-Rice fc Lada I (119911 ) have argued that 



the distribution of high velocity CO gas in this main outflow can only be explained by a 
superposition of two independent outflows with different orientations, one almost head-on 
and the other more inclined to produce the bipolar appearance. Despite the overwhelming 
emission from these two main outflows, two additional possible outflow s associated w i th the 
Mon R2 cluster have been discovered by CO emission line mappin g by iTafalla et al. Jl997h 
and a third possible outflow (a microjet from IRS 3) was found by Preibisch et al. ( 20021 ). 



The 12 CO(2-l) and 13 CO(2-l) velocity channel maps of ITafalla et all (119971 ) (their Figs. 
3 and 4) show a redshifted high velocity compo nent (Mon R2-N) about 80" to the north of 
the central cluster, coinciding with "Object 1" of ICohen &: Frogel I (119771 ). This object is the 
brightest of a group of embedded stars north of the ring-like reflection nebula, and is itself 
embedded in reflection nebulosity with a complex filamentary structure shown in Fig. 2. 
The brightest star i n this sub-group of the Mon R2 cluster is located at 6 h 7 m 45?9, -6°21'47". 
Also. ITafalla et alJ <jl997h discuss a strongly blueshifted 12 CO feature (-100",-40") from IRS 1 
that also coincides with an enhancement of CS emission, indicating the presence of a deeply 
embedded object. We refer to this object by the initials of the authors as "TBWW97 CO- 
100,-40". At its position (6 /l 7 m 39-2, -6°23'44"), we note a faint, slightly elongated infrared 
object on the S(l) image, surrounded by several small patches of nebulosity (Fig. 2 and 
Fig. 3, upper panels). 

The Spitzer 4.5 /mi image shows the same object morphology and confirms that the 
marginally extended 2.12 /im object is not an embedded star. The Spitzer 8.0 /j,m image 
shows extended emission in this region. The high-velocity CO gas and the shock-excited H 2 
emission indicate a spatially compact outflow, without a clear jet morphology. This object 
is shown near the western edge of Fig. 2 (labeled "CO"), and in the first panel of Fig. 3. 



The infrared source IRS 3 was noted by (IBeckwith et al.lll976l ) to be extended and was 
later resolved into a num ber of sources. The hi ghest spatial resolution image are the speckle 
interferometry results by lPreibisch et al.l (120021 ) that list point sources IRS 3 A - F, of which 
A, B, and C are associated with strong reflection nebulosity, and B is p ossibly associated with 
a microjet of about 0.5" length. The detailed CO velocity maps by iGiannakopoulou et al. 



-10- 



( 119971 ) show that IRS 3 is not associated with the big Mon R2 CO outflow, but that it is 
a separate source of very high-velocity CO outflow, even though a clear bipolar feature was 
not detected. 



3.3. The Spatial Distribution of H 2 Jets 



The newly discovered H 2 jets are overlayed on Fig. 1 as open circles with labels. 
Al so overlayed are t he lowest contours of the blue and redshifted features in the CO map 
of IWolf et all Jl990h The ir m ap was chosen over the higher spatial resolution maps of 
Meyers-Rice fc Lada I (Il99ll ) and iTafalla et al.l (119971 ) since it covers a larger area with higher 
sensitivity to faint extended features. The only sites of star formation lying clearly outside 
of the CO outflow lobes are the newly discovered HOD07 1 to the west of the main cluster, 
and the previously known young object GGD 11. The longest of the newly discovered H 2 
jets (HOD07 13) lies just south of the redshifted CO lobe. 

The CO outflows probably have a sh ell structure with a relatively empty cavity near the 
outflow axis (jMeyers-Rice &: Lada Ill99ll ). Since the outflowing material interacts turbulently 
with the molecular material of the ambient cloud, triggered star formation would be expected 
near the interface between the shell and the ambient cloud. In projection, the highest density 
of triggered star formation sites would be expected where we look tangentially along the shell 
wall, and a smaller number where we look onto the front and back sides of the outflow shells. 
The overall distribution of newly found H 2 jets outside of the central cluster roughly matches 
this expectation. Note in particular that the sma ll cluster as s ociate d with HH 866 is projected 
against the tip of the blueshifted CO contour of I Wolf et al.l (jl990l ). and that the group of H 2 
jets HOD07 2, 3, 6, 8, and 9 lie close to the eastern edge of the blueshifted CO contour, and 
behind two secondary blue-shifted CO emission maxima. These spatial coincidences strongly 
suggest that sites of recent star formation outside the main Mon R2 cluster may have been 
triggered by turbulent interaction of the massive main Mon R2 outflow with the ambient 
molecular material. 



In the following, we discuss the newly discovered H 2 jets individually on the basis of 
their morphology and, when available, by a comparison of the near-infrared 2.12 /mi S(l) 
image to longer wavelength Spitzer images. The objects are numbered in the sequence of 
increasing R.A., and labeled by the prefix HOD07. All coordinates are given in the J2000.0 
system. To allow an easy comparison of their sizes, the overview 2.12 /j,m images based on 
the UKIRT/WFCAM data are all shown on the same spatial scale. Similarly, the detailed 
views in the figures comparing WFCAM and Spitzer/IRAC results are all on the same, but 
finer, spatial scale, with the exception of Fig. 6, where this was not practical. 
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The small localized reflection nebulae that we found on our 2.12 /zm image, but did 
not classify as jets, are listed in Table 2 and are indicated by small squares in Figs. 1 and 
2. These localized reflection nebulae are concentrated to the south of the Mon R2 cluster, 
closer to the cluster than the H 2 jets found to the south and south west of the cluster. Since 
the reflection nebulae represent a more advanced state of star formation than the jets, this 
distribution might indicate a sequence of triggered star formation to the south of the main 
Mon R2 cluster, in the region where the blue and red shifted CO contours overlap. 

3.4. Discussion of Individual H 2 Jets 

3.4.1. HOD07 1 

The H 2 jet HOD07 1 has a very symmetric bipolar appearance, as is shown in Fig. 4. 
We note a high degree of similarity in the shape and relative intensity of corresponding 
knots in the two lobes of the jet. We interpret the two extended patches of nebulosity 3" 
east (labeled IE) and west (labeled 1W) of the apparent center of symmetry at coordinates 
Qh 7 m 10 s _ 6 o 26 ' 34 " ( mar ked by a circle in Fig. 4) as emission or scattered light from the two 
outflow cavities, and assume that the driving source is hidden from direct view and located 
near that position in the dark region between those two emission knots. From this position, 
the eastern jet lobe has a position angle of ^ 60°. Since HOD07 1 is not included in the 
publicly released Spitzer images of the Mon R2 region, this position of the central source 
could not be confirmed by longer wavelength data. 

Going symmetrically outwards from the assumed position of the driving source, there 
is a pair of faint emission knots at 6 ft 7 m 9-5, -6°26'44" (2W) and 6 ft 7 m ll-6, -6°26'22" (2E), 
respectively. The next pair of relatively faint emission knots is at 6 ft 7 m 9-0, -6°26'46" (3W) 
and 6 ft 7 m 12-4, -6°26'16" (3E). Further away from the central source we find two complex 
systems of bright emission knots, containing two (west) and three (east) individual knots, 
centered at 6 ft 7 m 8-2, -6°26'52" (4W) and 6 A 7 m 12.9, -6°26'11" (4E), respectively. Both of 
these emission regions show a bend of ~ 45° to the left in the direction of motion. To the 
west of the central source, at 6 h 7 m 7*5, -6°26'54" (5W) there is one additional very strong 
emission knot of approximately triangular shape with no detectable counterpart in the east. 
Finally, the most distant, approximately symmetric set of faint shock emission features is 
located at 6 A 7 m 6-3, -6°26'57" (6W) and 6 ft 7 m 16-7, -6°25'53" (6E) On the eastern side, this 
faint extended nebulosity coincides with the images of several stars, but there is no evidence 
to suggest that these stars are physically associated with the outflow. 

It is noteworthy that the bend in the otherwise straight bipolar jet coincides with the 
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strongest emission knots. Bent or S-sha p ed je ts as a result of disk pr ecession in a binary 



system were discussed by iTerquem et al.l ( 119991 ) and iBate et al. I (120001 ). We speculate that 



the main accretion and outflow event that produced these strong systems of emission knots 
was triggered by a close stellar encounter, most likely with a binary star component on 
an elliptical orbit, and that this encounter was also responsible for a change in the disk 
orientation and therefore the jet direction. 



3.4.2. HOD07 2 

The extended S(l) emission with shock-like morphology labeled HOD07 2 in Fig. 5 at 
6 h 7 m 30-5, -6°11'51" consists of two distinct knots of emission, suggesting that this is a short 
jet (P.A. fh 120°) originating from a star that itself is too obscured to be directly visible. 
The position of the more compact shock emission knot is the position given above and in 
Table 1. The Spitzer 8 /mi image (Fig. 6, bottom panel) shows extended emission centered 
on a point source at 6 ft- 7" 1 30-9, -6°11'48". While this point source indicates other ongoing 
star formation activity in this region, its position does not coincide with the jet HOD07 2 
and it therefore cannot be identified with the driving source of this jet. 



3.4.3. HOD07 3 

The bow shock HOD07 3 (Fig. 5) at 6 /l 7 m 31-6, -6 12'47" is probably driven by a source 
some distance away to the south east of the shock. Projecting back from the bow shock in this 
direction (P.A. « 147°) we find a string of faint, slightly extended knots of 2.12 fim emission. 
These faint emission features, indicated by arrows in Fig. 5, are located at 6 h 7 m 33^2, - 
6°13'03", at 6 h 7 m 34-l, -6°13'24", and at 6 h 7 m 39-3, -6°14'01". Overall, they appear to form a 
jet with multiple, faint internal shocks, ending in the bright bow shock. All these emission 
knots are confirmed by the Spitzer 4.5 /txm data, even though only the knot closest to the bow 
shock is included in Fig. 6, where we show the bow shock at 2.12 (xm, 4.5 /zm, and 8.0 /im. 
However, the position of the driving source of this jet could not be identified, neither on the 
WFCAM image nor the Spitzer images. We speculate that the jet emerges from deep within 
the molecular cloud and that only the terminating bow shock moving towards the observer 
is seen at relatively low extinction. Therefore, as indicated in Fig. 5, the position listed for 
this jet in Table 1 is that of the bow shock. 
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3.4-4- HOD07 4 

Two knots of H2 emission are visible to the west of the star at 6 /l 7 m 35*3, -6°20'00" in 
Fig. 7 and the top panel of Fig. 8. However, the extended nebulosity associated with the 
south-western of these knots is not connected with this star, strongly suggesting that this 
star may not be physically associated with the H 2 emission. A second, fainter system of H 2 
emission with a bow shock morphology is visible at 6 ft 7 m 35^9, -6°19'56", to the east-northeast 
of this star. The IRAC 3.6 /im (not shown here) and 4.5 /mi images (Fig. 8, top center) show 
faint curved nebulosity connecting the 2.12 /im features with a very faint star at position 
6 /l 7 m 37*6, -6°19'54" that is visible on the 2.12 /mi UKIRT image as well as on the two shorter 
wavelength IRAC images. There is no evidence for a counterjet, bent or straight, emerging 
from the position of that faint star. Unfortunately, the longer wavelength IRAC archival 
images at this position are not of useable quality. From the available data we conclude that 
the S(l) shock emission and the S-shaped feature in the IRAC 3.6 /xm and 4.5 /tm images 
represent the less obscured parts of a jet emerging from the star at 6 /l 7 m 37^6, -6°19'54". To 
explain the bending of the jet, it must be emitte d from a precessing d isk, similar to the case 



of IRAS 03256+3055 in NGC 1333 discussed bv lHodapp et all fl2005h . 



3.4.5. HOD 07 5 

The position given here for HOD07 5: 6 /l 7 m 38-7, -6°21'14", indicated by a circle in Fig. 7 
and Fig. 8, middle panel, is that of a faint extended emission feature near the apparent 
symmetry center of a system of shocks oriented at a position angle of fh 32°. This central 
object is visible on the UKIRT 2.12 /im image, but is much brighter on the Spitzer 4.5 /tm 
image, indicating heavy obscuration. Also, the 2.12 /im image of this central feature shows 
a marginally detected bifurcation, indicating a disk seen in absorption against the extended 
emission or a disk shadow. We tentatively conclude that the driving source is not directly 
visible at near-infrared wavelengths, but that emission and/or scattered light very close to 
the source has been detected, and that the collimating disk is seen nearly edge-on. From 
this driving source, the jet extends along a position angle of ~ 30°. The closest pair of shock 
emission knots is about 5" to the north and south of the driving source. A second pair of 
emission features ^20" to the north and south appears less symmetric, with the southern 
shock extending further away from the source. 
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3.4-6. HOD 07 6 

To the west of the center of Fig. 9 and in Fig. 10 (top panels), a faint string of four 
2.12 /iin S(l) emission knots (labeled 6) and one larger emission knot further to the south 
west at 6 /l 7 m 38*7 -6°11'59" are all roughly aligned at P. A. m 45°. The position of the central 
source cannot be determined with any certainty, and the Spitzer images do not provide 
additional information in this case. For lack of any better information, we list this jet by the 
coordinate of the approximate center of symmetry of the four northernmost emission knots 
located at 6 /l, 7 m 41^2, -6°11'14", implicitly assuming this as a plausible position of a central 
source if these knots are, as in many other cases, symmetric shock fronts in a bipolar jet. 

3.4.7. HOD07 7 

The position given for this object in Table 1 is that of a faint star at Q h 7 m Al -5, -6°21'28" 
roughly between two bright knots of S(l) emission in Fig. 7 and Fig. 8 (bottom panels) that 
define a jet axis along P.A. m 30°. Faint filamentary emission extends back from the shocks 
to the immediate vicinity of this star, strongly suggesting that it is the source of the jet. 
On the 2.12 /im image, the northern jet bends to the east and the southern jet to the west 
about 8" - 10" from the star, suggesting a strongly precessing jet. The two brightest S(l) 
emission knots are located near the bend in the jet (in the north) and slightly beyond it in 
the south. This jet morphology can be understood by assuming that the bent, contiguous 
jet was formed by a precessing driving source. The two emission knots lying outside of the 
contiguous jets could be the product of a single burst of higher outflow velocity. The fact 
that these higher velocity knots were emitted in the same projected direction as the strongest 
bend in the jet may indicate that the strong precession of the driving source's accretion disk 
that caused this bend may also have caused the higher accretion rate that in turn produced 
the higher outflow intensity and velocity found in the two brighter knots. 

The IRAC 3.6 /iin (not shown here) and 4.5 /im (Fig. 8, bottom center) images show 
the strongest, and spatially extended emission associated with the northern of the two S(l) 
emission regions. To the east of this object, extended emission is seen both on the S(l) image 
and on the IRAC 4.5 /iin and 8.0 /iin images, and this emission appears to be associated with 
a filament extending from the Mon R2 central cluster. Whether this filament is physically 
close to the driving source of the outflow or whether the outflow interacts with the filament 
cannot be conclusively determined. However, the different emission flux levels in the northern 
and southern lobe could be understood if the northern jet ran into denser material. 
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3.4-8. HOD 07 8 

On the 2.12 /j,m image in Fig. 9, this jet extends from an unresolved object at 6 /l 7 m 43*7, 
-6°10'46" to a knot at 6 A 7 m 44-4, -6°10'11", a position angle of 16° and distance of 36", 
corresponding to a projected length of 0.14 pc. The jet axis is outlined by extended 2.12 /jm 
emission. The only bright 8.0 /zm object (Fig. 10) associated with this jet is the southern 
emission knot. It is therefore likely that this 8.0 /zm object is the central star driving a jet 
where predominantly only one side is visible. The position of this bright mid-infrared source 
is listed as the position of the jet in Table 1 and shown in Fig. 9 and Fig. 10. On the 2.12 /zm 
image, emission extends about 2.5" to the south of the unresolved central source, suggesting 
that this is part of the jet's southern lobe. 

3.4.9. HOD 07 9 

In Fig. 9, a well developed bow shock is located to the west of a star at 6 A 7 m 44-0, 
^ll^". An emission feature behind the bow shock, near the center of the arc formed 
by the bow, could be a Mach disk associated with the bow shock. The driving source of 
this bow shock cannot be determined with certainty from the data available to us. We are 
discussing two possible identifications for the driving source: The most plausible candidate 
for the driving source of this bow shock is the star at 6 /l 7 m 44*0, -6°11'02" that is the second 
brightest 4.5 fim and 8.0 fxm source in Fig. 10 and that shows extended emission on the 
S(l) image indicating an illuminated outflow cavity, opening towards the west. We use this 
position at 6 A 7 m 44-0, -6°11'2" as the nominal position of this jet in Table 1. 

A much fainter, slightly extended object is visible on the UKIRT S(l) image at 6 7 m 43*6, 
-6°10'58" and is also detectable, partly blended with the object described just above, on the 
Spitzer 4.5 ^m image. The only reason to consider this object as an alternative candidate 
for the driving source of the bow shock is that it is more precisely located on the symmetry 
axis of the bow shock. This object, on the other hand, is also located on the axis of the jet 
source HOD07 8, and may just be another shock front of the southern part of that jet. The 
brightest source in Fig. 10 at any of the wavelengths discussed here is a cometary reflection 
nebula at 6 /l 7 m 43-4, -6°11'16" that we list as reflection nebula R5 in Table 2. 

3.4.IO. HOD07 10 

Just west of its center Fig. 2 shows 2.12 /im S(l) emission with the morphology of a jet 
ending in a bow shock at 6 /l 7 m 43-0, -6°24'20". This object is also shown in the bottom panels 
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of Fig. 3. The bow shock is strongly detected in the Spitzer IRAC 4.5 /tm image and is also 
indicated in the 8.0/xm image. It is not clear where the driving source for this jet is located. 
The jet can be traced to an extended feature located at 6 /l 7 m 44*l, -6°23'59" and detected in 
both the 2.12 /im and IRAC 4.5 /im images. It is embedded in the emission surrounding the 
central Mon R2 cluster. This may be just another knot in a jet emerging from deeply inside 
the molecular cloud, or this may indicate the position of the driving source. There is no 
enhanced flux in the 8.0 /im image at that position. Nevertheless, we adopt the position of 
this putative driving source as the formal coordinates of this jet called HOD07 10. A second 
shock front, located at 6 fe 7 m 43-6, -6°24 / ll" appears as the brightest of a series of internal 
shocks along the jet axis, which is oriented at a position angle of 220°. 



3.4.II. HOD07 11 

The region shown in Fig. 11 contains 2 separate systems of H 2 knots shown in more 
detail in Fig. 12. A well-defined bow shock located at 6 A 7 m 48-6, -6°19'03" and several fainter 
knots of 2.12 /im S(l) emission south-west from it suggest a jet with an axis oriented at P.A. 
~ 23°. The last S(l) emission knot at 6 fe 7 m 46-9, -6°20'6" coincides with a faint Spitzer 4.5 
/im object (Fig. 12, center panel), but is most likely just the 4.5 /im emission of the same 
emission knot, and not the driving source itself. The Spitzer 4.5 /im and 8.0 /im images in 
Fig. 12 show a stronger mid-infrared source at approx. 6 /l 7 m 46?8, -6°20'12" just south-west 
of the position of the last observed S(l) knot in this string. This object is not detected in 
our near-infrared S(l) image and therefore must be very deeply embedded. We tentatively 
identify this mid-infrared source as the driving source of the jet HOD07 11. 



3.4.12. HOD07 12 

The object at 6 A 7 m 49-8, -6°20'43" in Fig. 11, HOD07 12, has the morphology of a faint 
bipolar nebula seen nearly edge-on, suggesting that it could be the source of a jet. The 
Spitzer 4.5 /im image in Fig. 12 shows this object as an unresolved source, while the 8.0 /im 
image does not show the object. Two other knots of S(l) emission lie to the south-east 
at 6 A 7 m 51-6, -6°20'58" and at 6 ft 7 m 54-9, -6°21'57", respectively, suggesting a jet axis at a 
position angle ^152°. While there are a number of faint emission features to the north-east 
of the bipolar nebula, e.g., at 6 /l 7 m 48*8, -6°20'23", no clear association of these features with 
the jet source can be established. We tentatively identify the bipolar nebula as the driving 
source of HOD07 12 based on its morphology. A much brighter star, also associated with 
extended emission at 2.12 /im and in the Spitzer bands, and located at 6 /l 7 m 49*14, -6°20'33", 



-17- 



lies also on the jet axis and is surrounded by some 2.12 jum emission. It is the brightest 
source in the area on the Spitzer 8.0 /im image, indicating that it is a luminous, deeply 
embedded object. 

3.4.13. HOD07 13 

The nearly symmetric system of H 2 shocks labeled HOD07 13 and located south of the 
Mon R2 cluster is the largest jet found in our survey. The center of the outflow HOD07 13 
is clearly identifiable by an object on our S(l) image (Figs. 13 and 14) at position 6 h 7 m 57-4, 
-6°31'6". From this object, the northern jet extends along a position angle of ~ -13°. Unfor- 
tunately, the available 4.5 /im and 8.0 /im Spitzer images do not extend sufficiently far south 
to include this object. However, all 2.12 /im features of this object are confirmed by the 
Spitzer 3.6 /tm image (Fig. 13, right panel) that shows essentially the same outflow shock 
features. 

A magnified view of the central region of this outflow is shown in Fig. 14. Immediately 
north of the central object, a faint string of emission knots is visible stretching about 15" 
to the north. There is no counterpart to this knotty jet to the south of the central source, 
indicating the the jet is inclined with its northern part towards the observer, and that 
the southern counterjet is obscured by molecular material surrounding the central source. 
However, the IRAC 3.6 /im image shows a faint outline of a paraboloid outflow cavity wall 
to the south. 

More distant S(l) emission shock regions show a remarkable level of symmetry, since 
they are not affected by the obscuring material surrounding the central source. The first 
pair of emission regions are two complex regions, each showing two strands of emission 
knots, centered roughly at 6 fe 7 m 56-3, -6°29'39" (labeled IN) and 6 A 7 m 57-9, -6°32'23" (labeled 
IS). Further out, we find a symmetric pair of fainter shock fronts at 6 /l 7 m 55*7, -6 29'4" 
(labeled 2N) and 6 /l 7 m 58*3, -6°32'56" (labeled 2S). The shock fronts farthest removed from 
the central source show substructure that is again similar in both the northern and southern 
jet: A symmetric pair of shock fronts at 6 ft 7 m 54-l, -6°27'33" (labeled 3N) and 6 A 8 m 00 ? 2, 
-6°34'30" (labeled 3S). Finally, the leading shock fronts are found at 6 ft 7 m 53-8, -6°27'19 
(labeled 4N) and 6 /l 8 m 01-0, -6°34'43 (labeled 4S). The total projected length of this jet is 
w7.5', corresponding to 1.8 pc. 
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3.4-U- HOD07 14 

The compact, bipolar jet HOD07 14 (Fig. 15 and Fig. 16 top panel) is characterized by 
two bright shock fronts, at 6 /l 7 m 57-8, -6°25'26" and 6 fe 8 m 00-0, -6°25'40", respectively, that 
define a jet axis at position angle ~ 114°. The Spitzer 4.5 fim image (Fig. 16, top row, 
center) confirms the emission seen at 2.12 fim, but also shows additional emission, possibly 
indicative of a second jet parallel to HOD07 14, but more deeply embedded and therefore 
not visible in the near-infrared. The only indication of the position of the driving source of 
HOD07 14 is a flux minimum in the Spitzer 8.0 /im image, possibly a dense molecular clump 
seen in absorption against background emission. This argument would place the driving 
source at approximately 6 /l 7 m 58'^9, -6°25'31", the position listed in Table 1 and marked by a 
circle in Fig. 16. 



3.4.15. HOD07 15 

The UKIRT S(l) image (Fig. 15) shows a nearly symmetric, bent system of shock fronts, 
with two bright shocks both to the north and south of the assumed position of the driving 
source. In addition, on each side, one fainter, extended shock is visible further away from 
the center of symmetry. The Spitzer 4.5 //m image (Fig. 16, bottom row, center) shows two 
additional emission knots in the jet, closer to the center, while the Spitzer 8.0 /im image 
shows a flux minimum in between those two strongly obscured 4.5 /mi shocks. We identify 
this position, approximately 6 h 8 m 10^0, -6°24'47", as the position of the driving source of the 
jet HOD07 15 in Table 1 and Fig. 15. 



3.5. Other Bipolar and Cometary Nebulae 

The area near the Mon R2 cluster contains numerous small patches of nebulosity, often 
with bipolar or cometary morphology, that did not fit our criteria for identification as shocks 
or jets. We list those objects in Table 2 and indicate them by small square symbols in 
Figs. 1 and 2. The bipolar or cometary shape of many of these objects suggests that these 
are young, embedded stars still surrounded by disks and that they have just excavated an 
outflow cavity in the surrounding molecular material. These objects appear particularly 
numerous in the area about 2 to 3 arcminutes south of the main Mon R2 cluster, and just 
north of t he cluster, a s is dem onstrated in Figs. 1 and 2. This finding is consistent with the 



result by iGutermuth I (120051 ) that the Class I sources in Mon R2, identified by their J, H, 



and K, and Spitzer IRAC colors, are concentrated in a filamentary distribution to the south 
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of the cluster center, and to the north and north east of the cluster. While there is no strict 
relationship between reflection nebulae of bipolar or cometary morphology and SED Class 
I, the two criteria cover objects of similar evolutionary status at the trailing edge of their 
accretion phase. The main diff erence between the distribution of reflection nebulae and of 
the Class I sources identified by lGutermuth I (120051 ) is that we find fewer individual reflection 
nebulae in the region dominated by the large, ring-like reflection nebula and the filamentary 
emission that dominates the center of the cluster. We attribute this difference largely to the 
difficulty of identifying small reflection nebulae around individual stars by morphological 
criteria in this crowded central region of the Mon R2 cluster. 



By comparison, iGutermuth I (120051 ) found that the older Class II sources, identified by 
J, H, and K-band colors, are more uniformly distributed around the Mon R2 cluster center, 
indicating a dynamically more relaxed population of more developed young stars. 



OTHER INTERESTING OBJECTS IN THE FIELD 



4.1. HH 866 



The group of Herbig-Haro objects HH 866 found by lWang. Stecklum. fc Benning I (120051 ) 
in the L 1646 cloud was included in our field (Figs. 1 and 17). The HH objects ar e associated 
with the IRAS source 06046-0603. The region was identified by lCarpenter I (120001 ) as a region 
of enhanced s tar de nsity in 2MASS data, and was identified as a potential CO outflow source 
bv IXu et all (l200ll ). Our UKIRT/WFCAM 2.12 pm image shows a system of features that 
can be morphologically identified as likely H2 shocks. 

The main system of emission knots in the center of Fig. 17 is extended roughly in 
north-south direction, and is probably physically associated with the bow shock at 6 /l 7 m 07-5, 
-6°04'36". The position of the driving source of this jet cannot be determined with certainty. 
We tentatively list the center of symmetry of the brightest H2 emission knots at 6 h 7 m 07^9, 
-6°03'42" as the position of the driving source, and name this jet HH 866 Jet W in Table 1. 
The central parts of th is north- south jet are coinciding close ly with knots A, B, and C in the 
Ho; and [SB] images of IWang. Stecklum. fc Henning I (120051 ). To the east of this position is a 
large, rather poorly defined bow shock at 6 h 7 m 08-l, -6°03'36" that appears to be associated 
with more shock emission knots further east of it. A plausible identification of the driving 
source is a star embedd ed in nebulosity at 6 /t 7 m 09^8, - 6°03'4 4". The nebulosity around this 
star was also noted by IWang. Stecklum. fc Benning I (120051 ) and the star is labeled " 1" in 
their Ba image. Our image resolves this object into two stars, the eastern and fainter of 
which is associated with the nebulosity and is the likely outflow source. We list the position 
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of this star as HH 866 Jet E in Table 1. 



4.2. GGD 11 



The bipolar reflection nebula GGD 11 (iGyulbudaghian. Glushkov. fc Denisyuk Ill978l ) 



was included in our UKIRT/WFCAM image and in Fig. 1. We show a more detailed view of 
this object in Fig. 18, because our image is, to our knowledge, the best available near-infrared 
image of GGD 11. This object has the typical morphology of a bipolar nebula, with the 
eastern lobe being far brighter than the western lobe, due to inclination of the object. The 
brightness distribution is fairly smooth, in contrast to the knotty and filamentary appearance 
of typical H 2 jets. 



4.3. NGC 2182 

The NGC 2182 reflection nebula was included in our UKIRT/WFCAM image even 
though it is not included the image section shown in Fig. 1. It is shown in detail in Fig. 19. 
Different from the smooth appearance of this object on optical photographs, e. g. on the red 
plate of the DSS2, and on the K-band 2MASS image, the UKIRT 2.12 /mi S(l) image shows 
two filamentary regions of enhanced S(l) line emission to the west of the illuminating star. 
We cannot determine the excitation mechanism of these two filaments just from one image 
and both fluorescent excitation and shock excitation are possible. 



5. COMPARISON WITH OTHER STAR FORMING REGIONS 



In our 2.12 /im imaging survey of the Mon R2 molecular cloud for H 2 jets associated 
with forming stars, we have found a total of 15 new H2 jets in Mon R2 and two additional 
iets in the L 1646 (HH 866) region. Together with the two outflows det ected in CO 
by iTafalla et al.l (119971 ) and the IRS 3 micro jet found by iPreibisch et al.l (120021 ) this gives a 
minimum of 20 active jets and outflows in the larger area of Mon R2 (includ ing the L 164 6 
region), in addition to the huge, but probably inactive, CO outflow found by lLoren I (Il98ll ). 
The spatial distribution of the H 2 jets in Mon R2 along the projected edge of the large CO 
outflow, in particular the cluster of H 2 jets 12' north of the Mon R2 main cluster, suggests 
that this star-formation activity was triggered by the huge fossil CO outflow in Mon R2, one 
of the largest outflows known. 



By comparison, in a similar survey of Orion A, 



Stanke. McCaughrean. fc Zinnecker 
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( 120021 ) found 44 jets with a high degree of certainty, and 29 less certain cases. Considering 
that their survey area was larger than ours, and the Orion A cloud is more massive than the 
Mon R2 cloud, this comparison implies a similar specific star formation rate in both clouds. 
The level of outflow activity in Mon R2 i s also simila r to th at found in the two massive 
star-forming regions DR21 and W75, where iDavis et al.l (120071 ) found a combined number of 
approximately 50 H2 jet sources based on a similar UKIRT WFCAM. 

NGC 1333 is one of the mo st active star-forming regions in the solar vicinity and 
Bally. Devine. fc Reipurth I (119961 ) list 15 optical Herbig-Haro objects, H 2 jets, and CO out- 
flows in this region. This region is closer to the Sun than the other regions discussed before, 
and the discovery of jets is therefore easier. In the P erseus molecular cloud as a whole, of 
which NGC 1333 is a part, IWalawender et al.l (120051 ) counted a total of 141 optically (Ha 
and [SII]) detectable Herbig-Haro objects that, at the minimum, belong to 30 individual out - 
flows. The smaller Barnard 1 cloud (~ 1200 M^) was found by IWalawender et al. I (l2005bl ) 
to have 8 protostars driving outflows. Also, iBally et al.l (120061 ) found a total of at least 20 
outflows in the Chamaeleon I molecular cloud, mostly by an optical search supported by 
Spitzer data. Chamaeleon I is a nearby (165 pc), relatively low-mass (1000 M ) cloud. In 
light of the different techniques used to search for collimated jets from young stars in these 
different regions, we conclude that the number of such jets discovered in Mon R2 is similar 
to that found in the other massive star forming regions. 

While it may be argued that the presence of reflection nebulosity, fluorescently excited 
H 2 emission and PAH emission makes the discovery of shock-excited H 2 emission near the 
central cluster of Mon R2 difficult, this is also the region most intensely studied in the past, 
and only the IRS 3 microjet was found there. By contrast, most of our newly discovered 
jets are located in the periphery of the central cluster, and the longest collimated jets were 
found farthest away from the central cluster. The complete absence of large H 2 jets near 
the center of the Mon R2 cluster cannot be explained by observational selection effects alone 
and suggests that H 2 jets do not form in this environment. 



This is in agreement with the finding by lDavis et al.l (120071 ) in DR21/W75 who conclude 
that the H 2 jets in their sample are usually not associated with young molecular cores 
detectable at sub-mm wavelengths or with compact infrared clusters. They further conclude 
that clustering may inhibit disc accretion and the production of extensive flows. A similar 
conclusion can also be drawn from the exte nsive survey of the Orion A cloud for H 2 jets by 
Stanke. McCaughrean. fc Zinnecker I (120021 ). In Orion A, H 2 jets ar e predominantly found 
outside the area of the Trapezium cluster and the Orion nebula. Also, iHodapp I (119941 ) noted 
in their K' imaging survey of CO outflow sources that member objects of large young clusters 
are less likely to show jets, or localized reflection nebulosity than young objects in smaller 
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groupings. 

6. Conclusions 

Our UKIRT WFCAM imaging survey of the Mon R2 molecular cloud for H 2 jets as- 
sociated with forming stars was based on a 2.12 /im H 2 1-0 S(l) plus continuum image. 
Object identification was based on morphology and confirmation by archival Spitzer images. 
We have found a total of 17 new H 2 jets, including the two in L 1646, as well as 27 small 
reflection nebulae. 

We conclude that the number of H 2 jet sources in Mon R2 is similar to that found in other 
star forming regions of comparable size. We note that the H 2 jets appear to be concentrated 
in projection against the walls of the outflow cavity created by the large Mon R2 CO outflow. 
This is indicative of a scenario of sequential star formation in Mon R2 triggered by the large 
outflow. 

We confirm the finding in other high-mass star forming regions that H 2 jets tend to be 
found outside of the dense central cluster in such regions. It appears that the environment 
of a dense cluster containing some massive stars prevents the formation of large H 2 jets. 

The United Kingdom Infrared Telescope is operated by the Joint Astronomy Centre on 
behalf of the U.K. Particle Physics and Astronomy Research Council. 

This work is based in part on observations made with the Spitzer Space Telescope, which 
is operated by the Jet Propulsion Laboratory, California Institute of Technology under a 
contract with NASA. 

This publication makes use of data products from the Two Micron All Sky Survey, 
which is a joint project of the University of Massachusetts and the Infrared Processing and 
Analysis Center / California Institute of Technology, funded by the National Aeronautics 
and Space Administration and the National Science Foundation. 
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Fig. 1. — Image of the Mon R2 star forming region in the 1-0 S(l) emission line of H2 at 
2.12 /iin, obtained with the WFCAM at UKIRT. The circles indicate the positions of the 
newly found H 2 jets. The Herbig-Haro object HH 866 is visible in the north-west corner 
of the image. The small reflection nebula nea r the eastern edge of the image is GGD 11 
(jGyulbudaghian. Glushkov. &: Denisyuk I Il978l ). Newly found small reflection nebulae are 
indicated by squares. Superposed on this image are lowest contours of the blueshifted (-2 to 
6 km s — 1, solid lin e ) and redshifted ( 14 to 22 km s— 1, dashed line) CO emission from the 
map by I Wolf et all Jl99ph . 




Fig. 2. — A detailed view of the Mon R2 center region, based on the same UKIRT/WFCAM 
2.12 /im S(l ) line + continuum image as Fig. 1. The two CO outflows identified by 
Tafalla et al.l (119971 ). Mon R2 N and CO 100,-40, are indicated, as well as the newly found 
out flows HOD07 7 and 1 0. The positions of the three dominant infrared sources, IRS 1 
- 3 ( Beckwith et al.lll976l ). are indicated in the overexposed central region of the Mon R2 
cluster. The positions of many of the small reflection nebulae from Table 2 are indicated by 
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Fig. 3. — This figure shows detailed views of the two outflows in the Mon R2 center region 
in the 2.12 /im S(l) line + continuum, Spitzer IR AC 4.5 //m, and Sp itzer IRAC 8.0 /iin 
bands. The top panels show the outflow CO 100,-40 (ITafalla et al.lll997l ). the bottom panels 
the outflow HOD07 10. 
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Fig. 4. — This 2.12 fj,m S(l) line + continuum image shows outflow HOD07 I. The tentative 
outflow center (marked by a circle) is located at 6 /l 7 m 10-4, -6°26'34" in the flux minimum 
between the two patches of nebulosity that define the center of symmetry. The symmetric 
sets of shocks in the western and eastern jet are labeled. 
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Fig. 6.— This figure compares the 2.12 ^m, IRAC 4.5 fim and 8.0 /im images of two 
separate systems of H 2 knots: HOD07 2 at Q h 7 m 30-5, -6°11'51" is a small jet with several 
internal shocks and the position given is that of the most compact knot. The 8.0 /xm Spitzer 
image shows extended nebulosity and a point source that, however, cannot be identified as 
the central source of the jet. The bow shock HOD07 3 is located at Q h 7 m 31-Q, -6°12 / 46". 
The driving source of this bow shock could not be identified. 
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Fig. 7. — This 2.12 S(l) line + continuum image contains 3 separate systems of H 2 
knots: HOD07 4 at 6 A 7 m 35-3, -6°20'0", appears as a star with nebulous knots. As shown 
in Fig. 8, this position is most likely not the central source of this outflow. HOD07 5 at 
6 ft 7 m 38-7, -6°21'14", is the likely central star of this bipolar jet. HOD07 7 at 6 A 7 m 41-5, 
-6°21'28", is the likely position of the outflow center of that jet. 
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Fig. 8. — This figure contains the 2.12 /im S(l), 4.5 and 8.0 fim images of the three 
systems of H 2 knots found in the area of Fig. 7: The driving source of HOD07 4 is probably 
located at Q h 7 m 37-Q, -6°19'54", marked by a very faint star and the apex of the jet-like 
nebulosity of the Spitzer 4.5 /mi image. HOD07 5 lies at 6 /l 7 m 38-7, -6°21'14", identified by a 
faint source at the approximate center of symmetry of the bipolar jet. HOD07 7 at 6 fe 7 m 41*5, 
-6°21/28", is the likely position of the driving source of a bipolar jet. 
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Fig. 9. — This 2.12 /im S(l) line + continuum image contains 3 separate systems of H 2 
knots: HOD07 6 at 6 h 7 m 41*2, -6°11'14" is the center of symmetry of a string of emission 
knots. HOD07 8 at 6 ft 7 m 43-7, -6°10'46" is the likely center of a single-sided jet. HOD07 9 
at 6 A 7 m 44-0, -6° 11' 2" is a star associated with H 2 emission and the likely source of a bow 
shock. 
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Fig. 10. — This figure shows the S(l), 4.5 ^m, and 8.0 /xm images of three systems of H 2 
knots: In the top row, HOD07 6 at 6 ft 7 m 41-2, -6°11'14" is the center of symmetry of a string 
of emission knots. In the bottem row, HOD07 8 at 6 fe 7 m 43^7, -6°10'46" is the likely center 
of a single-sided jet. HOD07 9 at 6 h 7 m 44*0, -6°11'2" is a star associated with H2 emission 
and the likely source of a bow shock. 
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Fig. 11. — This figure is a 2.12 /zm S(l) line + continuum image and contains two separate 
systems of H 2 knots: The circle identifying HOD07 11 at 6 ft 7 m 46*9, -6°20'6" marks the 
putative position of the central source of a jet extending to the north-north-east (P.A. ^23°) 
and ending in a well-defined bow shock located at 6 ft 7 m 48?6, -6°19'3". Several fainter knots 
of S(l) emission lie between the center and the bow shock, completing the morphology of 
a typical jet. The object HOD07 12 at 6 fe 7 m 49-8, -6°20'43" has the morphology of a faint 
bipolar nebula. Two other knots of S(l) emission lie to the south-east at 6 h 7 m 51-6, -6°20'58" 
and at 6 /l 7 m 54-9, -6°21'57 // , respectively. While there £1X6 db number of faint features to the 
north-east of the bipolar nebula, no clear association with this jet source can be established. 
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Fig. 12. — This figure shows images in 2.12 /zm S(l), 4.5 yum, and 8.0 /im of the central 
region of the two outflow sources shown in Fig. 11: HOD07 11 at 6 A 7 m 46*8, -6°20'12" and 
HOD07 12 at 6 A 7 m 49-8, -6°20'43". 
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Fig. 13. — This figure shows the largest of the newly discovered jets: HOD07 13 with the 
central source located at 6 ft 7 m 57-4, -6°31'6" The jet has a projected length of 7.5', corre- 
sponding to 1.8 pc. The left panel shows the UKIRT / WFC AM 2.12 fim S(l) + continuum 
image, the right panel is the Spitzer/IRAC 3.6 ^m image. The area of HOD07 13 was not 
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Fig. 15. — This 2.12 /im S(l) line + continuum image contains two separate systems of H 2 
knots: The likely central source of HOD07 14 is located at 6 ft 7 m 59-2, -6°25'33". The central 
source of HOD07 15 is located near 6 A 8 m 10-0, -6°24'47". 
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Fig. 16. — This figure contains the 2.12 /zm, Spitzer 4.5 yum, and 8.0 /im images of two 
separate systems of H2 knots: The likely central source of HOD07 14 is located at 6 /l 7 m 59?2, 
-6°25 / 33". Based on the symmetry of the emission features, the central source of HOD07 15 
located near 6 A 8 m 10-0, -6°24'47". 
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Fig. 17. — This 2.12 /iin S(l) line + continuum image of the nebulosi ty in L 1646 associated 
with the optical HH object 866 jWang. Stecklum. fc Henning I boOfih at 6 h 7 m 7-8, -6°3'49". 
We identify two H 2 jets labeled Jet W and Jet E in Fig. 2 and listed in Table 1. 
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Fig. 18. — This 2.12 /iin S(l) line + continuum image shows the bipolar reflection nebula 
GGD 11. 
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Fig. 19. — This 2.12 fim S(l) line + continuum image shows the reflection nebula NGC 
2182, which was recorded in our large WFCAM image but not included in Fig. 1. In contrast 
to optical and infrared continuum images, this image shows strong, filamentary and knotty 
S(l) emission suggesting emission from the shock front of an expanding shell. 



-45 - 



Table 1. Outflows in Mon R2 



RA (J2000.0) 


Dec (J2000.0) 


Name 


6 h 7 m 07-9 


-6°03'42" 


HH 866 Jet W 


6 h 7 m 09-8 


_6°03'44" 


HH 866 Jet E 


6 h 7 m 10-4 


-6°26'34" 


HOD07 1 


6 h 7 m 30?5 


-6°11'51" 


HOD07 2 


6 h 7 m 31-6 


-6°12'47" 


HOD07 3 


6 h 7 m 37-6 


-6°19'54" 


HOD07 4 


6 h 7 m 38-7 


-6°21'14" 


HOD07 5 


6 h 7 m 39-2 


-6°23'44" 


TBWW97 CO -100,-40 


6 h 7 m 41-2 


-6°11'14" 


HOD07 6 


6 h 7 m 41 ? 5 


-6°21'28" 


HOD07 7 


6 h 7 m 43-7 


-6°10'46" 


HOD07 8 


6 h 7 m 44-0 


-6°11'02" 


HOD07 9 




-6°23'59" 


HOD07 10 


6 h 7 m 45-9 


-6°21'47" 


TBWW97 Mon R2-N 


6 h 7 m 46-8 


-6°20'12" 


HOD07 11 


6 h 7 m 49?8 


-6°20'43" 


HOD07 12 


6 h 7 m 57-4 


-6°31'06" 


HOD07 13 


6 h 7 m 58-9 


-6°25'31" 


HOD07 14 


6 h 8 m 10-0 


-6°24'47" 


HOD07 15 
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Table 2. Reflection Nebulae in Mon R2 



RA (J2000.0) 


Dec (J2000.0) 


Name 


Morphology 


6 h 7 m ll-0 


-6°26'51" 


HOD07 Rl 


reflection neb. 


6 h 7 m 36 f 6 


-6°15'07" 


HOD07 R2 


bipolar neb. 


gh 7 m 41 ? 3 


-6°20'02" 


HOD07 R3 


comctary neb. 


6 h 7 m 43?3 


-6°18'27" 


HOD07 R4 


reflection neb. 


6 h 7 m 43 f 4 


-6°11'16" 


HOD07 R5 


cometary neb. 


6 h 7 m 43-9 


-6°19'30" 


HOD07 R6 


cometary neb. jet ? 


6 h 7 m 44 f 3 


-6°25'40" 


HOD07 R7 


cometary neb. 


6' i 7 m 44-8 


-6°20'31" 


HOD07 R8 


reflection neb. 


6 h 7 m 45-5 


-6°25'44" 


HOD07 R9 


reflection neb. 


6 fe 7 m 45-5 


-6°26'08" 


HOD07 RIO 


reflection neb. 


6 h 7 m 45-9 


-6°25'29" 


HOD07 Rll 


comctary neb. 




-6°25'31" 


HOD07 R12 


reflection nebula 


6 h 7 m 46-8 


-6°25'07" 


HOD07 R13 


cometary neb. 


6 fe 7 m 46 fg 


-6°24'17" 


HOD07 R14 


bipolar neb. 


6 h 7 m 48-l 


-6°24'05" 


HOD07 R15 


reflection neb. 


6' i 7 m 48-2 


-6°25'07" 


HOD07 R16 


reflection neb. 


6 h 7 m 48?6 


-6°25'22" 


HOD07 R17 


bipolar neb. 


6 h 7 m 48?8 


-6°25'55" 


HOD07 R18 


cometary neb. 


6 h 7 m 50-7 


-6°25'53" 


HOD07 R19 


reflection neb. 


6' i 7 m 50-7 


-6°26'10" 


HOD07 R20 


reflection neb. 


6 h 7 m 51 ? 


-6°23'00" 


HOD07 R21 


bipolar neb. 


6' i 7 m 51-l 


-6°26'10" 


HOD07 R22 


reflection neb. 


6 h 7 m 51-3 


-6°19'47" 


HOD07 R23 


bipolar neb. 


6 h 7 m 51 f 5 


-6°26'23" 


HOD07 R24 


reflection neb. 


6' i 7 m 51-7 


-6°23'18" 


HOD07 R25 


reflection neb. 


6 fe 7 m 51 f 8 


-6°25'40" 


HOD07 R26 


reflection neb. 


Qh 7 m 54 ? 7 


-6°26'36" 


HOD07 R27 


reflection neb. 



